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Abstract 

For the family of 0(n) invariant nonlinear a-models we con- 
sider boundary conditions that are periodic up to an 0(n) ro- 
tation. The helicity modulus is related to the change in free 
energy under variations of the corresponding angle. It defines 
a nonperturbative finite volume running coupling similar to the 
Schrodinger functional for QCD. For the two-dimensional 0(3)- 
model we investigate this quantity by analytical and numeri- 
cal techniques. We establish its universal continuum relation to 
the finite volume massgap coupling at all scales and coupling 
strengths. 
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1 Introduction 



The property of asymptotic freedom is a decisive feature of QCD as well 
as of a large class of two dimensional nonabelian spin models like the O(n) 
cr-models for n > 2. Although it is based only on proofs in perturbation 
theory (to all orders), the following structural properties of these models 
are widely accepted and assumed hereQ. The continuum limit is reached 
at vanishing bare coupling and a mass scale emerges in the renormalized 
continuum theory by dimensional transmutation. Many features associated 
with much higher energies or short distance can be related to each other 
by renormalized perturbation theory. Usually one singles out one suitable 
high energy quantity as renormalized coupling and constructs expansions for 
other observables in its powers. Other phenomena around the fundamental 
scale, like the spectrum, are not accessible to this strategy and are often 
investigated by numerical simulation. These opposite "sectors" are really 
features of one and the same theory and it is hence both interesting and 
possible to relate them, that is, to compute renormalized coupling constants 
at high energy in terms of the low energy scale. This has been the programme 
of the ALPHA collaboration in recent years. An efficient method has been 
developed first for the 0(3) a-model || followed by quenched QCD which is 
reviewed in [|]]. Many present activities are related to the goal of including 
dynamical quarks. 

To relate the perturbative sector with low energy physics very dissimilar 
continuum scales have to be accommodated on a lattice with a spacing that 
is small compared to all other scales. In a direct approach this either calls for 
unmanageably large lattices or one has to compromise with the attempted 
limits like the continuum extrapolation. The ALPHA strategy overcomes 
this problem by a finite size scaling technique. A running coupling constant 
g 2 (L) is constructed in the continuum limit, which at large L can be related to 
spectral scales and which at small L can be used as an expansion parameter 
and thus gives access to the perturbative sector. Step by step one computes 
g 2 (2L) in terms of g 2 (L) by continuum extrapolation. Since the system size 
L is used as the physical scale, L/a is the only large scale ratio, where a is 
the lattice spacing. The choice of g is not unique and a number of practical 
criteria were taken into account. 

2 See ref. [|| for a diverging point of view 
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For the O(n) model the finite volume mass gap was used in || 

f = ^-m(L)L, (1.1) 
n — 1 

where m(L) is the gap of the transfer matrix referring to spatial periodic 
boundary conditions with periodicity L. For L — > oo, m(L) saturates to the 
infinite volume massgap, which we identify with the dynamically generated 
scale. At small L the mass gap becomes perturbative || and can be used as an 
expansion parameter. For gauge theories a similarly convenient though less 
obvious quantity was defined via the Schrodinger functional |], Q . The basic 
mechanism is to study the response of the free energy under the variation of 
an angle that enters into non-trivial boundary conditions. Again the system 
size is the only scale beside the cutoff, and for the quenched theory a s in the 
high energy limit has been computed to a satisfactory precision. 

In this article we define and investigate an alternative coupling gr{L) for 
the O(n) cr-model in two dimensions, which is closely related to the helic- 
ity modulus. We compare its properties with the massgap coupling g(L). 
They can be analytically related in the continuum for both small and large 
coupling. Both expansions are checked and the crossover range is controlled 
by high precision numerical simulations. In the next section we define the 
helicity modulus and relate it to g at strong coupling. In section 3 gx is 
properly normalized and its weak coupling behavior is explored up to two 
loops with details given in appendix A. Section 4 summarizes our numerical 
work and gives conclusions. This work is based on the diploma thesis |l] 
of H. Molke. The introduction of gy goes back to earlier attempts |7|, to 
investigate renormalization by finite size techniques. 

2 Helicity modulus and transfer matrix 

In this section we introduce the helicity modulus. For earlier discussions of 
similar quantities see Ref. and further references found there. 

We consider the O(n) cr-model with its standard nearest neighbor lattice 
action ^ 

S — — s(x)s(x + afi), (2.1) 

where s(x) is the unit length spin field and /t, fi — 0, 1, are unit vectors along 
the axes of a square lattice with spacing a. We take T as the size in the time 
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or O-direction and L in the other direction. In space direction we impose 
strictly periodic boundary conditions, while in the time direction we demand 
periodicity up to a planar SO(ra) rotation in spin space by an angle a, 

s(x + TO) = e^(aKij)s(x). (2.2) 

Here generates rotations in the ij plane, 

( Kij )ki = SikSjl ~ SilSjk ' ( 2 '^ 

and we assume i ^ j fixed to an arbitrary pair of values until further notice. 
Integration of all spins with the 0(n) invariant measure gives the partition 
function 

Z a = J Dsexp(-S). (2.4) 

Ratios Z ai /Z a2 depend on differences in free energy for different boundary 
conditions and are expected to be universal. We now define the helicity 
modulus T by 



T = _ J_ d2Z « 



(2.5) 

a=0 



Z da 2 

From this definition it is rather easy to establish a connection between T 
and the transfer matrix T. For the fully periodic case Z we have 

Z = ti[T T ' a }. (2.6) 

The extra twist by an angle a corresponds to inserting a rotation operator 
under the trace in the Hilbert space where T operates, 

Z a = tr [T T/a exp(ia/Cij)]- (2.7) 

If we realize states as wavefunctions ip(cr) on spatial one dimensional spin 
fields <j, this induced operator is given by 

(^exp(ialCij)ipj(a) — ^(exp(— aKij)a^. (2.8) 

The operator T possesses real positive eigenvalues A > Ai > A 2 > . . . 
with Ao corresponding to the nondegenerate ground state. These eigenvalues 
depend on L and a. We define the finite volume massgap m(L) as 

exp(-m(L)a) = ^. (2.9) 
Ao 



Due to the 0(n) invariance T and /C^- commute and we simultaneously diag- 
onalize the generator that appears in (|2.7| ). Hence for each eigenstate there 
is a value = 0, 1, 2, — For n = 3 these are just the "magnetic" quan- 
tum numbers of eigenstates in all possible integer isospin multiplets in the 
spectrum. The a-dependent partition function is now given by 

Z a = £(A fc ) T/a exp(za/i fc ), (2.10) 

fc>0 

and for the helicity modulus we obtain 

T = ^E(^) T/ >L (2.H) 
z o fc > 

For large mT we expect the ground state and the lowest excitations above 
it to dominate. While the ground state is O(n) invariant, /x = 0, we expect 
an n-fold degenerate vector multiplet of one particle states as the next exci- 
tations above it. On it the generators are represented in the form ( |2.3| ) and 
have eigenvalues 1,-1 and n — 2 times 0. Asymptotically this implies 

X T ~°°2]Texp(-£ p T). (2.12) 

p 

The sum here is over single particle momenta p = 2irj/L with —L/2a < j < 
L/2a, and E p are the associated energies. Contributions of higher states are 
assumed to be exponentially suppressed. We now take the continuum limit 
at fixed large values for mL and mT. The one-particle spectrum is expected 
to become relativistic, E p = y/ m 2 + p' 2 and we find that T is given in terms 
of the renormalized coupling (|1.1| ) introduced in @ by 



]=-oo 

with 



T ~ 2 exp -PV7 2 + (27TJ) 2 • (2.13) 



n-l_ 2 



7 = — (2-14) 
This formula holds in the continuum limit for large g 2 and fixed aspect ratio 

p=~. (2.15) 
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Figure 1: Relative discrepancies between the strong coupling formula (|2.13|) 
and approximations to it. 



For not extremely large g 2 it only takes the few lowest j-values to carry out 
the sum in ( [2.13| ) numerically to machine precision. By Poisson resummation 
one may derive the asymptotic form 



T 



27 



ir 1 (p7) + 0(exp(-^p 2 + l 7 )) 

^exp(- P 7)(l + 0( 7 - 1 )), 
np 



(2.16) 
(2.17) 



where K\ is a modified Bessel function. In Fig. 1 we show how these asymp- 



totic forms approach the sum (|2.13| ) for the case p = 1. Also shown is the 
j = term alone, which underestimates the sum by at most 5 % for 7 < 3.5. 
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3 Helicity modulus in perturbation theory 

3.1 Preparation and leading order 

By changing the integration variables in 



we arrive at 



where the action 



s(x) -> exp ( -a^Kij ) s(x) ( 3.1 } 



Z a = J Dsexp(-S B ) (3.2) 

S B = — \ s(x)B fl s(x + apt) (3.3) 
9o xfi 



has acquired a constant SO(n) gauge field 

B = exp U^K^ , B 1 = l, (3.4) 

and s(x) has become strictly periodic. The helicity modulus is now expressed 
by an expectation value in the periodic ensemble 

with B^ = 1 after taking all derivatives. So far the angle a has referred to 
rotations with one particular generator Kij. At this point we average over all 
planes i < j and split T into two O(n) invariant contributions T = Ti — T 2 , 

M^}^£? <sWs(i+o6)> (3 - 6) 



and 



'dS B V\ 2 a 4 



In these expressions the volume V = TL and the currents 

tfAx) = s{x)K l0 ^s{x) (3.8) 
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with the discrete derivative 

A^s(x) = - [s(x + afi) — s(x)] (3.9) 

have been introduced. While Ti is a nearest neighbor correlation essentially 
corresponding to the internal energy, T2 is a kind of susceptibility with corre- 
lations over all separations. In perturbation theory contributions to Y2 start 
at the one loop level and we find to leading order in g 

r = -^ + o(i). (3.10) 

A correctly normalized and nonperturbatively defined renormalized L-dependent 
coupling constant can now be defined as 

9 2 r = ~L (3.H) 
np T 

quite in the spirit of the Schrodinger functional. Its relation to g from ref. 
is 

2 _ J 9 2 + 0(^ 4 ) for f - 

9r 1 ^v^^y exp ^ps 2 ) for 9 2 -> 00. 

As g 2 is proportional to L at large volume or strong coupling, we find expo- 
nential growth for g\. This is also expected for the Schrodinger functional 
coupling in gauge theory flO |. The origin in both cases is the exponentially 
small sensitivity of the free energy to boundary effects in a physically large 
volume. 



3.2 Results of a two loop calculation 

For simplicity we confined ourselves to the case T = L, p = 1 for our per- 
turbative and our numerical calculations. Details on the perturbative eval- 
uation of T are given in appendix A. Nontrivial coefficients were evaluated 
for sequences of lattices of finite L/a and then fitted to the expected asymp- 
totic form. The extrapolation was carried out as described in appendix D of 



ref. [11| with lattices up to L/a = 100. The cost in CPU time was negligi- 
ble in this two dimensional case. As for the Schrodinger functional, it was 
advantageous to compute some of the two loop diagrams in position space 
rather than momentum space. 
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The nearest neighbor correlation in Ti is rather simple to compute \T2 
to the required order. As a result we find 

Ti = ~ - ^ ~ + O(g 4 , a 2 ). (3.13) 

n (?q 2n 16n 

The current-current correlation in T2 is more involved and has the structure 

T 2 = T( 0) + T^ 2 + O(^ 4 ). (3.14) 

By combining the results for and from the appendix with the 
expansion of g 2 on the lattice || |l6j we find in the continuum limit the 
renormalized perturbative series 

g\{L) = f(sL) + d x g\sL) + d 2 f(sL) + 0(f(sL)) (3.15) 

with 

d^s) = (n - 2) 0.16350689821 - — Ins (3.16) 

2n 

d 2 (s) = rfi(s) 2 + {n- 2) [0.00315826256 - (n - 2) 0.007733893180] 

' Ins. (3.17) 



(2k) 



Note that a free relative factor s between the scales at which g^ and g are 
taken has been introduced here for later use. 



4 Numerical results 

Our main goal in this section is to establish the nonperturbative relation 
between g 2 {L) and g\(L) in the continuum limit of the 0(3) model for ar- 
bitrary values of these couplings. Our strategy is to first construct series of 
values L/a and f3 = 1/ g 2 which correspond to fixed values of g. For precisely 
the same series we measure g\ and extrapolate these values to L/a — > 00. 
For the first part of the task we include but extend as necessary the data 
from p. These runs were carried out with precisely the same code as de- 
scribed there. In particular, we took advantage of free boundary conditions 
in the time direction to extract the massgap, and the reweighting technique 
allowed for a post-run fine-tuning of f3 to match a desired value of g 2 . The 
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Table 1: Simulation results for g 2 and g^. 



simulation of T is rather conventional on an L x L torus. We employed the 
single cluster algorithm [13| and measured the observables given in eqs. ( |3.6| ) 
and (|3.7|) . Both kinds of results are collected in Table |l|. 

For each of the eight series at fixed g 2 we have pairs of values with errors 
S of g 2 and (Sr) of g\. For the extrapolation in a/L we combine them into 
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Figure 2: Continuum extrapolation of g\ for g 2 
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(4.1) 



The required slope can be estimated with sufficient accuracy from the weak 
and strong coupling behavior. We then extrapolate by fitting a function A + 
B(a/L) 2 to the g\ values with these errors. A typical case is shown in Fig. 2 
where the coarsest lattice was not included in the fit. Recently there has been 
some debate in the literature [14, Q whether the use of an asymptotic cutoff 
dependence proportional to a 2 , based on perturbation theory, is justified. 
Our data are compatible with this form, and we assume it for continuum 
extrapolations in this work. As a variation, a linear cutoff dependence, as 
suggested for other quantities in |14| , could however also be fitted to the data 
in Fig. 2 and would lead to an exptrapolation to 4.377(7), which differs not 
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Table 2: Nonperturbative relation between g 2 and g\. 



dramatically but clearly significantly. To distinguish numerically between 
these and yet other behaviors would require much higher precision and is of 
interest for future study. All our extrapolated continuum results are collected 
in Table |^. They are plotted in Fig. 3 and the close-up in Fig. 4, where 
one sees the data branch off form the 2- loop curve ( [3.151) . After a very 



narrow transition region they approach the strong coupling behavior which 
we constructed by performing the sum in ( 2.13|) . 



To get a feeling for the lattice artefacts associated with our two couplings 
individually we estimated their step scaling function (SSF) for a pair of values 
corresponding to similar L, g 2 ~ 1.06 and g\ ~ 1.29. The SSF — the focus 
of interest in || - gives g 2 (2L) as a function of g 2 (L). Pairs of simulations 
at identical (3 and sizes L/a and 2L/a yield 

Y,{u ) a/L) = g 2 {2L)\ . (4.2) 

E is then extrapolated to the continuum limit. A completely analogous 
quantity Ex is defined in terms of g\. In Fig. 5 the continuum approaches 
of E and Ex are shown. Lattice artefacts amount to a few percent at say 
L/a = 8 with St deviating significantly more from its continuum limit than 
E. This trend is expected, since T receives contributions form several low- 
lying eigenstates of the transfer matrix (c.f. ( |2.11| )) while g is constructed in 
terms of the massgap only. 

Instead of relating the two couplings at one scale L one may also consider 
the connection between g\{L) and g 2 (sL). This has already been incorpo- 
rated in the perturbative formulas ( |3.15| ). One may hope that an appropriate 
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Figure 3: Coupling g\ versus g 2 with asymptotic expansions. 



choice of s improves the accuracy of the approximation as was for instance 



found in [10]. A somewhat natural value is s = 2.7936 (n = 3), for which cfi(s) 
in ( |3.16 ) vanishes and which coincides with the ratio of the A-parameters as- 
sociated with the short distance behaviors of the two couplings. To compare 
such an expansion with non-perturbative results we would have to simulate 
series of lattices of size L and sL at the same bare couplings and take the 
continuum limit. While this is not possible, we gained good control over the 
step scaling function for g 2 (L) in ||. We used its four loop approximation 
16| , [17| P1 and fitted the remainder with an effective five loop coefficient to 



evolve from g 2 (L) to g 2 (sL). In this way we found however no value s ^ 1 
which significantly improved the series for g\(L). 

To conclude, we have investigated two nonperturbatively defined cou- 
pling constants for the 0(3) nonlinear cr-model with exponentially different 

3 We evaluated b 4 = 0.0040 at n = 3 for eq.(3.47) in [jT|- 
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Figure 4: Close-up of g\/g 2 in the transition region. 

low energy behaviors. Analytical relations, valid in the continuum limit, 
are available for both weak and strong coupling. Precise numerical simu- 
lations covered the intermediate range and matched with both asymptotic 
expansions. As expected, the helicity modulus coupling shows larger lattice 
artefacts than the finite volume massgap. 

Acknowledgement We would like to thank Erhard Seiler for pointing out 
to us that the strong coupling asymptotics in an earlier version of this paper 
was incorrect. 

A Two loop expansion of T 

In this appendix we use lattice units with a = 1 and take T = L. For 
the perturbative expansion of the a-model on a finite lattice we have to 
fix the global 0(n) invariance by the Fadeev- Popov technique ||15|| . For 0(n) 
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Figure 5: Continuum approach of step scaling functions. 



invariant integrands I(s) it amounts to the replacement I(s) — > J(s)/(s)F(s) 
which does not change the value of the integral. For the noninvariant function 
/ we take 

f(s) = 6(M 1 )5(M 2 )...6(M n _ 1 ) (A.l) 

where M = J2 X s ( x ) i s the total magnetization, which is hence forced to 
point in the n-direction. The compensating Fadeev-Popov factor is in this 
case given by 

F(s) = (Ml™" 1 (A.2) 

up to an irrelevant overall constant factor. The spins are parameterized by 
ann-1 component real field <fr(x) 



The resulting contributions are gathered in an effective action 



(A.3) 



eflf. 



2 t— 1 



(A,0) 2 + 1(A mX ) 2 
9o 



[n 
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It is understood to be expanded in g , 

oo 

s cS . = J2g 2 k s {k \ (a.5) 

k=0 

and the function / still has to be included in the path integral. It leads to 
the omission of the zero momentum mode in the propagator of the field 
and makes perturbative coefficients now well defined. The last two terms in 
S c s, correspond to the O(n) invariant measure and to F. We shall only need 
the leading terms 

S (0) = ^E(A,0) 2 (A.6) 

= ^E(A/) 2 -^(l-(n-l^)D 2 (A.7) 
The term defines the propagator 

I , e ip{x-y) 

(M*)My))o = s k i G(x - y ) = 6 kl -J2 — ^— ( A - 8 ) 

' p p 

where 1 < k, I < n — 1 and the primed sum is over the appropriate lattice 
momenta excluding p = (0, 0) and we have introduced p fl = 2 sin(p /i /2). 
To compute Ti we set 

s(x)s(x + a6) = 1 - g 2 Q E^ - g 4 E^ + 0(g 6 Q ) (A.9) 

E (0) = ^(Ao0) 2 (A.10) 

E^ = i(A o 2 ) 2 (A.ll) 

o 

and find 

Ti = — 2 (l - g 2 (E^)o ~ 9t [<£ (1) >o - (^ (0) 5 (1) )o]) (A.12) 

with the connected correlation in the last bracket. These contributions eval- 
uate to 

(£<»))„ = ^(l-^) (A.13) 

<£<»}„ -W>$ = !^i(i- ; l) a - ("- 1 )("- 2 ) ^.( A .i4) 

15 



The contribution T 2 has been expanded in ( |3.14 ). We introduce 
1 



, 4 EE %(*)m = V {0) + &BF> + 0(g*) 

V % i<j xy 



H <P) 



V 



E E ^ A o<; 



k<l \ x 



with the symmetric derivative 



and find 



T, 



n(n — 1) 
Numerical values are 



= o + A) - 00 - A)] 



(H(% = („_!)(„ -2)^--^) 
(# (1) >o = (n-l)[(n-2)S 1 + J B 2 ] 

(ffWfi^S = (n-l)(n- 2) [(^ - 1^(^ + 1^-1(1 

n - 2 , 

+ 2 A 3 

V 



1 



In the above expressions the following L-dependent constants were 
duced, 



At 
A 2 
A 3 

B 2 



Iv'— 



V 



(02)2 

Eu sin 2 ^ 



-^G(a:) 2 A M A / ,G'(.T) 



J B 1 -2^G(x)[A M G(x)] i 



A.15) 

A.16) 
A.17) 

A.18) 
A.19) 

A.20) 
A.21) 



A.22) 
intro- 

A.23) 

A.24) 

A.25) 
A.26) 
A.27) 



16 



Evaluated as x- or p-sums as they stand, only O(V) terms have to be summed. 
It is now straightforward to express the coefficients k\ , k2 in 



9t = — = fl'o + hgt + hgt + O(g 8 ) 



(A.28) 



in terms of the above constants. We evaluated them numerically for L = 
8 . . . 100 and determined the asymptotic behavior as explained in the ap- 
pendix of Ref . [ IT ] . The result is 



h 

k 2 -k\ 



(n-2) 
In -2) 



ML) 

2tt 

ML) 

(2vr) 2 



0.12165689529 
0.02514054821 



+ 



n — 1 



+ 0(1/L 2 ) (A.29) 



-On - 2) 2 0.00773389318 + — + 0(1/ L 2 ) 

96 



(A.30) 



Errors are beyond the digits given here, and the coefficients of ln(L)/L 2 and 
1/L 2 corrections are also known but not listed here. They are of the same 
size as the constants appearing here. To obtain the last fraction we set 
B2 = 1/48 + 0(1/L 2 ) which we observed to very high accuracy. 

The massgap coupling (1.1) has been computed perturbatively up to two 
loops in |3| and to three loops in JIB] , |T7| . From the last reference we extract 



9 2 = 9o+ migo + m 2 gl + O(^) 



(A.31) 



with 



mi 



m 2 — m 1 



(n-2) 
(n-2) 



ln(L) + ln(V2/7r)+7 



2tt 

ML) 

(2vr) 2 



2tt 

+ 0.021982285645 



+ i + 0(l/L 2 ) (A.32) 
+ ^ + 0(l/L 2 )(A.33) 



The exact fraction was again found to numerical precision. 
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